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Abstract 

We consider the quantum mechanical notion of the geometrical (Berry) phase in 
SU(2) gauge theory, both in the continuum and on the lattice. It is shown that in 
the coherent state basis eigenvalues of the Wilson loop operator naturally decompose 
into the geometrical and dynamical phase factors. Moreover, for each Wilson loop 
there is a unique choice of U(l) gauge rotations which do not change the value of the 
Berry phase. Determining this U(l) locally in terms of infinitesimal Wilson loops 
we define monopole-like defects and study their properties in numerical simulations 
on the lattice. The construction is gauge dependent, as is common for all known 
definitions of monopoles. We argue that for physical applications the use of the 
Lorenz gauge is most appropriate. And, indeed, the constructed monopoles have the 
correct continuum limit in this gauge. Physical consequences are briefly discussed. 



Introduction 



In his seminal paper ||T|, Berry found that a quantum state acquires an extra phase fac- 
tor when the Hamiltonian of the system depends on external parameters and the state 
is adiabatically transported along a closed path in the parameter space. The additional 
phase factor depends only on this path and reflects the geometry of the space of parame- 
ters. The mathematical construction behind the geometrical Berry phase was discovered 
by Simon 0, who showed that the adiabatic evolution of the quantum system naturally 
defines a line bundle structure over the parameter space. The geometrical phase is the 
holonomy of the U(l) connexion and cannot be removed by the global redefinition of the 
wave function when the bundle is non trivial. Moreover, the adiabatic approximation 
can be avoided and in fact any evolution of a quantum system naturally gives rise to the 
geometrical phase factor (for a detailed discussion of the non adiabatic Berry phase see 
and references therein). The Berry phase is again the holonomy of the U(l) connexion, 
but the base space of the bundle in question is now the space of all physical states. A 
unified consideration of the adiabatic and non adiabatic geometrical phases was presented 
in [§,0, where it was shown that these two, a priori different phase factors are actually 
related by the classification theorem for vector bundles. Note that the geometrical phase 
in quantum mechanics is ultimately related to the construction of the Abelian monopole 
since this is the only topologically non trivial object which arises when the structure group 
is U(l). 

In this paper we apply the notion of the Berry phase in the SU(2) pure gauge theory^ 
and, in particular, in the lattice formulation of the SU(2) gluodynamics. It is known that 
the geometrical phase construction is useful in chiral gauge theories (see, e.g., 0). Here 
we consider the Wilson loop operator in pure gluodynamics from the quantum mechanical 
point of view. Indeed, the P-exponent of the gauge potential is defined as a solution of 
the first order differential equation, which is very similar to the Schrodinger equation. 
Therefore the Wilson loop may be interpreted as the evolution operator of a quantum 
mechanical system, where the role of the time-dependent Hamiltonian is played by the 
gauge potential. Since the Wilson loop is an operator in the representation space of the 
SU(2) group, the space of the physical states is naturally described by the generalized 
coherent states 0, which parameterize the coset space SU{2)/U{1). We show that in the 
coherent state basis eigenvalues of the Wilson loop naturally decompose into dynamical 
and geometrical phase factors. For an infinitesimal Wilson loop we obtain the geometrical 
phase in terms of the continuum gauge potentials and argue that it reflects the non 
triviality of the Hopf bundle SU{2) SU{2)/U{1) = S\ 

Our considerations allow to define monopole-like defects in the SU(2) gluodynamics, 
which are different from the Abelian monopoles considered so far (see for a review). 
Indeed, the usual definition of Abelian monopoles implies a particular partial gauge fixing 
which leaves a U(l) subgroup unfixed and the monopoles are defined with respect to this 
remaining U(l). While in the standard treatment the U(l) is universally defined for all 

^ The SU(2) gauge group is taken for simplicity, the generalization to the SU(N) group is straightfor- 
ward. 
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the monopoles, in our construction the relevant U(l) subgroup is determined locally for 
each infinitesimal Wilson loop or, in the lattice formulation, for each plaquette. Therefore, 
the monopoles constructed are strictly speaking non-Abelian objects. The existence of 
the U{1) group inherent for each plaquette was in fact noticed in Ref. [|T^ and is central 
for the construction developed here. 

In particular, the consideration of the geometrical phase gives a new method to 
calculate the trace of the Wilson loop, which is periodic function of the phase angle, 
^TrW = cos(y9. We propose a method to calculate the angle ip without the restriction 
ip E [— vr; tt] and introduce the natural decomposition: 

(f = If mod 2tt + 2TTk , k E Z . (1) 

The Eq. (|1|), when applied to the elementary plaquette, gives a lattice analog of the Dirac 
strings associated with monopoles we have constructed. Therefore it becomes possible to 
consider the non-perturbative dynamics of these defects on the lattice. 

An evident problem, which can be understood already in terms of the quantum me- 
chanical analogy, is that the monopoles constructed in this way are not SU(2) gauge 
invariant. Following Ref. we argue that the most appropriate gauge is the Lorenz 



gauge, in which gauge dependent singularities of the gauge fields are maximally sup- 
pressed. By numerical calculations we show that the monopoles we have constructed are 
physical objects. In particular, the density of these monopoles is shown to scale correctly 
towards the continuum limit. We show that the asymmetry of the monopole currents 
behaves as an order parameter of the deconfinement phase transition. Finally, it is shown 
that these monopoles are locally correlated with the non-Abelian action density, like the 



monopoles in the Maximal Abelian gauge, Ref. ||12|| . 

The outline of the paper is as follows. In Section 1 we briefiy review some well known 
facts about the adiabatic and non adiabatic Berry phases in quantum mechanics. In Sec- 
tion 2 we consider the Wilson loop operator in the SU(2) gluodynamics. In Section 3 the 
SU(2) lattice gauge theory is considered. In particular, we present an explicit calculation 
of the Berry phase for an arbitrary Wilson loop on the lattice. In Section 4 we discuss the 
results of the numerical simulations. Our conclusions are summarized in the last section. 



1. Berry Phase in Quantum Mechanics 

Following the original paper of Berry consider a hermitian Hamiltonian H{t) = 
H{\i(t)) which is time-dependent through a set of parameters Aj(t), i = l,...,n. We 
are interested in evolution of the quantum mechanical system \ilj{t)) governed by the 
Schrodinger equation: 

tdtim) = Hixit)) \ m) (2) 

on the time interval [0;T]. Suppose, that for any fixed t G [0; T] the spectrum of H is 
non-degenerate: 

HiX) I n(A) ) = E„(A) I niX) ) ^ (3) 
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and that the time evolution of the parameters X{t) is slow, such that there are no in- 
duced transitions between energy levels. Then the adiabatic approximation applies and 
if the quantum mechanical system at t = is described by an eigenstate of Hamiltonian 
I -0(0) ) = I n{X{0)) ), the whole evolution reduces to a phase factor: 

\m) = e^^(*) |n(A(t))). (4) 

Moreover, the phase factor ip acquired by the state {tplt)) as a result of the adiabatic 
evolution consists of two parts [Q: 

t 

= e*7.W exp{-z j E„(r)dr} \n{\{t))) , (5) 



where S = — J En is the dynamical phase while the additional phase factor 7„ is known 
as the adiabatic Berry phase. The latter can be found by substituting (|) into the 
Schrodinger equation (^: 

dan = i{n{\)\^\n{\))dtX. (6) 
oX 

Consider now the parameters Aj(t) which are T-periodic, Aj(0) = Aj(T), and span a 
closed curve C in the parameter space. Then in the adiabatic approximation the initial 
and final state vectors differ only by a phase 

I V^(T) ) = e^^(^) I m ) . (7) 

A state vector which obeys Eq. (|^) is called the cyclic state. The Berry phase 'jniT) 
acquired by the cyclic state | ip ) may be represented in terms of the "gauge potential": 

A = ^{ni\)\^\ni\)), (8) 
7n(T) = [a = [ F , F = dA, (9) 



Sc 



where Sc is an arbitrary surface in the parameter space bounded by the contour C. In 
the form the geometrical nature of the Berry phase is almost evident. For this 

reason the Berry phase factor is also called geometrical phase. Indeed, the rich mathe- 
matical structure behind (|^, ^ was discovered by Simon and further elaborated in a 



number of papers p, ||, ||, [jT3|-[|T6[. In particular, it was shown in ||^ that the adiabatic- 
ity requirement is in fact unnecessary. Moreover, the precise form of the Hamiltonian 
is not important either. The geometrical phase factor depends only on the path C in 
the space of the physical states followed by the vector | ip ) and is naturally associated 
with any quantum evolution. For self-consistency we briefly review the essentials of the 
construction. 
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Consider unitary evolution of a quantum mechanical system described by a state vector 
\ilj{t)), where {ip) is an element of an + 1-dimensional complex vector space C^^^ 
with finite A^, for simplicity. In terms of the complex coordinates {zq, Zf^} the state 
vector is | ipit) ) = {zQ(t), ...,ziy{t)}. Since the unitary evolution preserves the norm, the 
normalization condition {ip \ = 1 defines a 2iV+ 1-dimensional sphere S*^^"*"^ G C^"*"^ on 
which the evolution takes place. The physical states are not given by normalized vectors 
m since | ip ) and {ip') define the same physical state if they differ by a phase 

\ip) = e^°'\ip')- Therefore, the set of physical states is the A^-dimensional projective 
space which is a Kahlerian manifold: 

CpN ^ 52A^+1/ f/(l) . (10) 

The quantum mechanical evolution is given by the Schrodinger equation (|^), where 
the Hamiltonian depends implicitly on time through the set of the parameters \i{t), 
Aj(0) = Xi(T). We are interested in a cyclic vector lipit)) which returns to the same 
physical state | ip{T) ) = e*"^*^"^^ | ipiO) ) after evolving along a closed path C G CP^. The 
phase (p(T) is the total phase acquired by the cyclic vector. Moreover, we introduce a 
single- valued vector \ip{t)) which differs from \ip{t))hj a. phase and satisfies the condition: 

|^(r)) = 1^(0)). (11) 

Note that Eq. (|TlD does not define the | ip ) uniquely: for a given \ip), vector e*^ \ip) is 
also single-valued provided that 9(T) — 9{0) = 2nn. 

In terms of the single- valued vector | 'ip ) the cyclic state | '^Z' ) is represented as: 

|^(t)) = e^^W \m) (12) 

and the Schrodinger equation gives: 

T 

^(T) = 5 + 7 = - j{i;\H\^) + t j{^\d\^). (13) 

c 

The first term in Eq. (|1^) is naturally identified with the non adiabatic dynamical phase 
and depends explicitly on the detailed structure of the Hamiltonian. One can show P 
that in the adiabatic limit it reduces to the expression 6 = — j En mentioned above. Note 
that the dynamical phase may also be calculated as 6 = — J {ip \H\ ip) since the phase 
difference between | tp ) and | ip ) drops out in the matrix element. The second term is a 
non adiabatic generalization of the geometrical Berry phase (^ and it depends only on the 
closed path C G CP^ spanned hj \-ip) during its cyclic evolution. Indeed, the geometrical 
phase 7 ' calculated by means of any other single- valued vector \ip' ) = e^^ \ip) (see the 
note after Eq. (pH])) differs from 7, Eq. (|13]): 

^' = i j{iP \ e-'^ de'^ \ ^) = -f- j dO = - 27m (14) 
c c 
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and therefore the difference 7 ' — 7 is inessential. 

In order to exphcitly evaluate the geometrical phase we note that single- valued vectors 
I '4) ) are parameterized by the homogeneous coordinates Wi on CP^ , \ip) = {wi, wn}, 
and the scalar product {tp' {ip) is given by: 

IV^) = {wi,...,wn}, I^') = {w[, ...,w'^} , 

(15) 

{^'\^) = exp{K{w',w) - lK{w',w') - lK{w,w)} , 

where w denotes complex conjugation and K{w\ w) is the Kahler potential on CP^ . One 
obtains [0, 0: 



7 = — / Im 



dK(w,w) , 
■aw 



dw 

c c 



1 r 1 

2 l + WiWi 



The appearance of a non trivial geometrical phase during a cyclic evolution of the 
quantum mechanical system is guaranteed then by the following topological arguments 
1^, |T^. Let variables Aj parameterize a manifold A. We have shown that the Hamiltonian 
H{\) defines the mapping / : A ^ CP^ which maps any curve in A to a curve in CP . 
bmce 7r„ (CP^) = 0, n ^ 2 and 7r2(CP^) = Z the mapping / might not be homotopic 
to zero provided that 7r2(A) 7^ 0. On the other hand, the line bundle (|10[) is non trivial 

and therefore the mapping 

g2N+i _^ CP^ ^ A defines a non trivial U(l) bundle over 
A. In this case the global definition of the phase of the state vector is impossible over 
parameter space and this fact guarantees the appearance of the Berry phase. Note that 
the geometrical phase factor when written in the form: 

7 = t f{i'\d\i') = I A = If, (17) 



Sc 



where Sq is an arbitrary surface which bounds the contour C G A, defines the winding 



number density j-P d\^{'4>\ d\.\il)) d'^X'^ of the mapping A 3 S"^ — > CP^ 



smce 



(d('i/' I) A (d| ^/^ )) is the first Chern class of the bundle (p^. In particular, the integral of 
^P over S*^ G A gives the degree of the mapping A 9 5^ CP^. This implies in turn 
that the integral ( p!7[ ) calculated for two different surfaces Sc and iS^ may differ only by 
27m and therefore the phase factor e^"' is surface independent. The non triviality of the 
mapping A 9 5*^ ^ CP^, together with a line bundle structure (p!OD, is essentially the 
only reason why the monopole potential repeatedly appears in studies of the Berry phase. 

Note that both dynamical and geometrical phase factors calculated in unitary trans- 
formed state basis are generally different. Indeed, going to unitary transformed states 
I ^Z'' ) = U\il)) is equivalent to going to the Hamiltonian H'{X) = U^H{X)U — iU^dtU. 
The corresponding total phase is then: 



T _ r _ _ ■ 

ip' = - J{iP\U+HU\ip) + i J {i>\U+dU\ij) + {^\d\ip) 
c '- 

= ~J{i,'\H\^P') + ^J{i,'\d\i,'). 
c 



:i8) 
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It is clear that both dynamical and geometrical phases calculated with H' differ in general 
from the phases evaluated with H. Nevertheless, there exists a class of unitary transfor- 
mations U (t) for which the Berry phases are the same. Namely, if for a unitary U (t) we 
have C = C, then 7 ' = 7 since the geometrical phase depends only on the path C G CP^. 
It is easy to see that only the U(l) transformations, U = e*" possess this property. 

2. Geometrical Phase and Wilson Loop 

There is a class of Hamiltonians which are of particular interest and are constructed as 
follows. The space of state vectors carries a unitary irreducible representation of a com- 
pact semisimple Lie group G. The Hamiltonian depends on a set of parameters Xi{t) and 
for every t is an element of the Lie algebra of G. Then the evolution operator for the 
Schrodinger equation (^, 

V{t) = Texp{-i [ H} (19) 
Jo 

is given by a path [0; T] ^ G in the group space and for any t belongs to the representation 
of G (in the Eq. (|l^) T denotes the time-ordering). Therefore, if at t = we start with 
an arbitrary state \ ip{0)), then the state vector at the time t, \ip(t)) = V(t)\4'{0)) is 
a generalized coherent state [§. It is important that the phase of the state \ip{t)) is 
naturally fixed with respect to a specific choice of the coherent state basis \z). Indeed, 
an arbitrary state vector may be represented as | ij{t) ) = e*^W| z{t) ) = e^^(*)[/(2(t))| ), 
where the complex variables {zi,...,Zn} parameterize the coset space G/H, H is the 
Cartan subgroup of G, U{z) G G/H and |0) denotes an arbitrary fixed vector in the 
representation space, which is usually taken to be highest or lowest weight vector. Thus 
the space of all physical states is nothing else but the coset space G/H. The total phase 
factor acquired by the cyclic state during the evolution (|^) may be calculated by means 
of Eq. dg) with = \z) 

The evolution operator (|19D is of particular importance in gauge theories. Indeed, 
consider a Wilson loop W{T) on the contour parameterized by coordinates x^{t), t G 
[0;r]: 

W{T) = Pexp{z [ A{t)dt}, A{t) = A^^^xit)) x^it) T% (20) 
Jo 

where T*^ are the generators of the gauge group G in the representation considered and 
the dot denotes differentiation with respect to t. The P-exponent is defined as a solution 
of the first order differential equations: 

(tdt + A)\i;) = 0, \m) = W{t)\m) (21) 

and therefore the Wilson loop is just the evolution operator ([19|) with a time-dependent 
Hamiltonian H = —A(t). Therefore, the discussion above suggests that phases of eigen- 
vectors of the Wilson loop may be naturally decomposed into dynamical and geometrical 
parts. Note that such a decomposition cannot be gauge invariant. Indeed, we have shown 
in the previous section that the dynamical and geometrical phases calculated for unitary 
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equivalent states vectors are generally different. On the other hand, the transformation 
H' = U^HU — iU^dtU is exactly the same as the gauge transformation. 

In the simplest case of G = SU (2) and the Wilson loop in the fundamental represen- 
tation, the procedure is as follows. The space G/H is a two-dimensional sphere which is 
parameterized by a complex coordinate z. The coherent states are build over the highest 
weight vector | ), cr^| ) = | ), a"' are the Pauli matrices: 



1 



(22) 



The action of the group element g G SU (2) in the coherent state basi^ is 



a j3 
—(3 a 



0(^,0 = arg[/5C + a] 



9\C) 



The Wilson loop operator, Eq. (^0]), has two eigenstates \z±): 



W{t) 



-(3{t) a{t) 



W{T)\z^) 



(23) 
(24) 

(25) 



It is sufficient to consider a particular eigenvector, e.g. | -z+ ), which tends to | ) when 
/? — s> (the Z- approaches infinity in the limit /5 — >■ 0). The cyclic state | C(^) ) a-^id the 
single valued vector | z{t) ) are constructed as follows: 



at) ) = Wit) I z+ ) 



z{t)) 



z{t) 



a{t) z+ - j3{t) 



(3{t) z+ + a{t) 

Therefore, the total phase factor f^iT) acquired by the cyclic state | C) is 

viT) = arg[/3(T)^+ + a(T)] = 5 + 7, 



(26) 



(27) 



5 



{z\A\z)= {C\A\C) 



7 



1^ Im 


r9/f(...) 1 




oz 



(28) 



where K{(,z) = ln[l + (z] is the Kahler potential on G/H. Similar to the quan- 
tum mechanical example, the dynamical phase may be calculated with either the single- 
valued vector or with the cyclic state. The trace of the Wilson loop operator is evidently 
iTriy(T) = cosifiT). 

^ Note that our notations are slightly different (although equivalent) to that of Ref. . 
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3. Lattice Implementation. 



Eqs. are directly applicable in the lattice gauge theories (LGT). Consider therefore 

the Wilson loop operator as it appears on the lattice: 



W 



9n 9n-i ■■■ 9i 



9i 



-Pi Oti 



Its eigenstate is determined by 



W\z+) 



—I 



Ima^ - sign(Ima„)A/(Ima^)2 + |/3^|2 



(29) 



(30) 



where the eigenvector which obeys the condition — >■ when — >■ has been se- 
lected. The calculation of the total phase is now straightforward. Indeed, according to 
Eqs. (E31E^ the cyclic | Cfc ) and single- valued | ) vectors are given by: 



k 



9k I Cfc-1 ) 
9k I Zk-l ) 

= 1 N, 



\Ck) 

zo = Co 



(31) 



Zn 



■tip 



z+. 



Therefore, the total phase factor ip of the Wilson loop is calculated as: 

N 

V = ^"Pk- 

k=l 



(32) 
reduces 



By construction, the total phase (p satisfies ^TrW = cos (p. Moreover, Eq. ([ 
to the well known Abelian expression when all Qi are diagonal. 

The decomposition of the total phase (|3^) into the dynamical and geometrical parts is 
more subtle. The point is that such a decomposition is not possible without a particular 
parameterization of the Wilson loop W and of the matrices gk in terms of a continuous 
time variable. We choose a natural parameterization: 



r 9iit) 



< t < 1 



W{t) = { gk{t)9k~i...9i, k-l<t<k 



(33) 



, 9N{t)...9i 



N -l<t<N 



where the time-dependent matrices 9k{t) are defined as follows: if gk = expjzAfc} then 
9kif) = exp{iAfc(t — k + 1)}. With the parameterization (|33D , the Wilson loop W(t) 
describes a continuous path in the SU(2) group manifold and the cyclic state | C(^) ) = 
W(t) \ Zj^ ) satisfies the Schrodinger equation idt\ ( ) = ~Ak\ (), k — 1 < t < k, k = 1, N. 
Since gk(t)Ak = A^gkit), we obtain from Eq. (03) the dynamical phase: 



N 



N 

s = ^(^fc-i 2fc_i) = y^(Cfc_i \Ak\Ck-i) 



(34) 



k=l 



k=l 
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where | ) are defined by Eq. ( plj) and the matrix elements ( z |A| 2;), A = yl"cr"/2 are 
given by 



[z\A\z) ^ 



^A^{1 - \z\^) + Re{zA-) 



A- = A^-iA^. (35) 



Moreover, the Berry phase associated with the Wilson loop (^) is 'j = ip — 6. 

Consider the Wilson loop which is the ordered product of the link matrices g^^^ taken 
along the boundary of the elementary plaquette situated at the point x and directed in 
the {/i, u} plane: 

W = 9t,u9t+u,|^9x+^,,u9x,,^ ■ (36) 
Let the area of the plaquette be a'^'^, then the dynamical phase S^^^uCr'^'^ is: 

^x,fiu C"^ ( ^+ I ^x,fi "I" 9x,fi-^x+fi,u9x,fi 9x,fj.9x+fj.,u^x+i/,fi9x+fi,u9x,fi 

~ 9t,fi9t+fi,u9x+u,tiAx,u9t+u,fi9x+fi,u9x,tJ. I ^+ ) o'^'^ ■ (37) 

Note that the state | z^ ) depends on the plaquette considered. For an infinitesimal Wilson 
loop Eq. (^) can be expanded in powers of a^'^. In the leading order the dynamical phase 
is given by: 

Sx,^uCt'''' = {z^\F^,-t[A^,A,]\z^)a^'' = { z+\D[^A,]\ z+) a^''' , (38) 

where F^^^, and are the continuum field strength tensor and covariant derivative, re- 
spectively: 

F,. = d[,A,] - t[A,, A,] = ^a^ {d[,A:^ + e^'''A\Al) , B^ = d,- i\A^, ] . (39) 

It can be checked directly that the total phase factor is ^x^iiv^^^^ = ( -2+ -2+ ) cr'^'^- 
Therefore, the Berry phase for the infinitesimal Wilson loop is given by: 

Ix^^.^y^'' = t{z+\[A„A,]\z+)a^''. (40) 

As expected from the quantum mechanical example above, the value of the Berry phase 
measures to what extent the gauge fields are non-Abelian. Indeed, for pure Abelian 
potentials = the Berry phase vanishes and the total phase of the infinitesimal 

Wilson loop is given by the dynamical phase alone: 

'Px,iJ.u = 8x,^u = ^dif,A^] . (41) 

As is mentioned above, the decomposition of the total phase of the Wilson loop into the 
dynamical and geometrical parts, Eqs. (|2^, ^D|), is gauge dependent. Nevertheless, 



it is possible to find a subgroup of the gauge group which does not change the Berry 
phase. Indeed, the defining property of these gauge transformations is the same as in 
the quantum mechanics, namely, they should not change the path spanned by the state 
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vector I z{t) ) in G/H = SU(2)/U{1). It is clear that the subgroup in question is a U(l) 
subgroup of SU(2). Since 1/2TtW = cosy?, the quantity (p mod 2n is gauge invariant 
and the effect of the U(l) gauge transformations which do not change the Berry phase, is 
to shift the dynamical or, equivalently, total phase by 27ik, k E Z. It is amusing to note 
that with k ^ these are indeed "large" U(l) gauge transformations familiar from the 
example of the U(l) compact electrodynamics. 

Moreover, when the gauge is fixed the considerations of the Section 1 are applicable. 
Namely, any closed path C in the physical space M is naturally mapped into a closed path 
in G/H, the relevant mapping being determined by the Wilson loop operator, calculated 
on C. The degree of the mapping / : 5 — G/H, where iS is a closed two-dimensional 
surface in M, is given by: 

deg[/] = ^1 1 = ^ I F, (42) 



27r J ' 27r „ 

S G/H 

where F is the first Chern class of the non trivial Hopf bundle G ^ G / H . Note that when 
deg[/] 7^ the Berry phase 7 contains a string-like singularity (Dirac string) somewhere 
on the surface S. As noted above, the position of the singularity may be arbitrary shifted 
by U(l) gauge transformations. 

Let the surface S be parameterized by coordinates x(a), a = {cr°} = {a^,a^}. Using 
Eqs. ( pot it is straightforward to show that: 

7 = - ^ y e'^'^ n-{a) ^(5) Al{x) dV^^ dV^, = e'^^d^x^d^x, , (43) 
5 s 



where the world-sheet vector field n°-(a), = 1 was introduced in Ref. IITO 



n [a 



(dV • F'') [(dV ■ F')^] , (dV ■ F") = dV^^F;, . (44) 

Note, however, that n"(o-), Eq. (|4^), is build on the F^^,, not = \e^u\pFxf, as in 

Ref. pi. 



4. Numerical Results. 

In this section we show that the formalism presented above may be used in lattice sim- 
ulations and present the numerical results which concern the total phase, Eq. (|3^) . The 
total phase of the Wilson loop is of particular importance since it naturally includes the 
contribution of both Abelian-like and non-Abelian fields. As for the lattice implementa- 
tion of the dynamical and Berry phases, we have not considered it yet although it seems 
straightforward. 

The trace of the Wilson loop, calculated over the boundary of an elementary plaquette 
p, is gauge invariant and is proportional to cosy^p. Therefore, the total phase may be 
represented in terms of the gauge invariant quantity (pp mod 27r and an integer number 
*p'- 

(Pp = (pp mod 27T + 27ik*p , (45) 
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where *p is the plaquette dual to p. By analogy with the corresponding equality in 
the compact U(l) LGT the integer k*p counts the number of Dirac strings piercing the 
plaquette p. A set of numbers k*p is an integer valued 2-form k G *C2{Z), which belongs 
to oriented plaquettes on the dual lattice. As usual, we define the monopoles as end points 
of the Dirac strings (^51) : 

= J e*Ci(Z). (46) 

The monopole current j^,^ is an integer valued 1-form on the dual lattice. While all the 
definitions above are unique in terms of the potentials (or link matrices), the potentials 
themselves are gauge dependent. As a result, the values of A;*p are gauge dependent as well 
and at first sight the definition (^) is devoid of any physical meaning. Clearly, the issue 
of the gauge dependence is crucial for physical applications of the procedure developed 
above. 

Therefore we introduce at this point a physically motivated gauge following the logic 
similar to that of the paper [|11|]. The observation is that in the continuum limit both 
the Dirac strings and monopoles correspond to singular gauge potentials A. It is easy 
to imagine, therefore, that by going to arbitrary large A, so to say inflated by the gauge 
transformations, one can generate an arbitrary number of spurious strings and monopoles. 
On the other hand, by minimizing potentials one may hope to squeeze the number of the 
topological defects to its minimum and these topological defects may have physical sig- 
nificance. In order to quantify what is understood by minimizing the potentials, consider 
the Lorenz gauge which is defined by the requirement that the functional 

R = ^(l-l/2Tr(7.,J (47) 

X,fJ, 

is minimal on the gauge orbit. In the naive continuum limit Eq. (^) reduces to -R = 
1/4 J(A^)^. Thus, the Lorenz gauge is singled out since in this gauge the link matrices 
are as close to unity as possible and the gauge dependent singularities are suppressed as 
much as possible. As we shall see in a moment, this heuristic justification of the gauge 
choice can be checked a posteriori. 

We have performed numerical experiments with monopoles defined by Eq. (HBI) in the 



Lorenz gauge. The calculations were done on the 12^ and 4 x 12^ lattices with periodic 
boundary conditions using 30 well equilibrated and statistically independent configura- 
tions. We used the local over-relaxation algorithm to fix the Lorenz gauge. The gauge 
was considered fixed when at each site the gauge transformation matrix fi^., which locally 
maximizes Eq. (0), satisfy 1 — l/2Trf2a, < 10~^. In order to circumvent the Gribov 
copies problem, each thermalized configuration was randomly gauge transformed to five 
gauge equivalent configurations and the Lorenz gauge was fixed again on each gauge copy. 



The measurements were performed on the configuration for which the functional (|47| ) is 
minimal. 



The density of monopoles, Eq. (46), is given by: 

1 

W 
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where V is the lattice volume. We have measured the density p in the Lorenz gauge as 
a function of the bare coupling (3 = The logarithm of p versus f3 on the symmetric 

12'* lattice is shown on the Fig. la. The solid curve on the figure is the renormalization 
group prediction: 



P 



const. ■ /3^53/i2i exp|_^^| ^ (49) 

plotted in the logarithmic scale. As is clear from the plot, the density p perfectly scales 
towards the continuum limit. Thus the monopoles (^Gf ) are physical objects and this 
justifies the choice of the Lorenz gauge, motivated above. Note that the finite physical 
density of the monopoles p6|) in the continuum gluodynamics implies in particular, that 
even in the Lorenz gauge the gauge fields are rather singular. The relevance of the singular 



fields in QCD is discussed in Ref. |19 



It is interesting to consider also the behavior of the density (^81) across the deconfine- 
ment phase transition. On the Fig. lb we plot the logarithm of p on the 4 x 12'^ lattice, the 
solid curve is the scaling law (^91) . The monopole density jumps at the critical coupling 
and remains non zero in the deconfinement phase. Note that after the phase transition p 
also seems to scale correctly towards the continuum limit, although this question needs 
further investigations. 

Another interesting quantity related to the monopole dynamics is the asymmetry of 
the monopole currents PUf: 



x,fi=l,2,3 X 

which is known to be the order parameter of the deconfinement phase transition when 
the monopoles are defined in the Maximal Abelian gauge. Indeed, the dual supercon- 
ductor confinement scenario suggests that the Abelian monopoles are condensed in the 
confinement phase, while they are almost static at high temperatures. Therefore, the 
asymmetry of Abelian monopole currents is unity at zero temperature and vanishes in 
the deconfinement phase with rising temperature. Since the definition of the monopoles 
( ^61) refers to the Lorenz gauge, the behavior of the asymmetry (^) across the phase 
transition is worth to be considered anew. 

The Fig. 2 represents the asymmetry A versus (3 on the 12^ (diamonds) and 4 x 12^ 
(squares) lattices. As expected, the asymmetry measured on the symmetric lattice is 
unity within the numerical errors. On the asymmetric lattice the asymmetry shows a 
rapid jump at the critical coupling, being substantially smaller after the phase transition. 
Thus, in the deconfinement phase the monopoles defined in the Lorenz gauge, Eq. (|i6|), 
are mostly static. Note that the asymmetry needs not to be unity at /? < jScrit on the 
asymmetric lattice since the Lorenz invariance is explicitly broken. 

In order to clarify the physical relevance of the monopoles constructed consider the 
correlation of the local SU(2) action density and the monopole currents [0. In particular, 
we define the relative excess of the action density associated with the monopole as follows: 

Sm — S 

V = g , (51) 
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where S denotes the average vacuum action density S =< 1 — l/2Trf/p >. The quantity 
Sm is the average value of 1 — 1/2 Tr ?7p calculated on the plaquettes, which belong to the 
3-dimensional cubes C^^f^ dual to the monopole currents j^^^: 



<l E (l-^Trf/,)>. (52) 



The average is implied over all cubes C^,^ for which j^^^ ^ 0. For static monopoles only 
the magnetic part of the SU(2) action density contributes to Sm- We have measured 
the quantity t] on the symmetric 12"^ lattice, the result is plotted on the Fig. 3. It is 
clearly seen that the monopoles in the Lorenz gauge are indeed locally correlated with 
non-Abelian SU(2) action density. 

Such an excess of the action associated with monopoles is quite a common phe- 
nomenon. Indeed, monopoles are correlated with instantons which realistically rep- 



resent non-perturbative background in QCD |22|. A possible general connection between 
monopoles and background fields was recently pointed out in Ref. |jl^. Moreover, the 
excess of the action associated with monopoles ( P6| ) is numerically close to the action 
excess around the Abelian monopoles in the Maximal Abelian gauge [|12 . 



The latter observation rises the question about the local correlation between these 
objects. A priori, it is by no means evident that the monopoles defined by the Eq. (^6]) and 
the Abelian monopoles in the Maximal Abelian gauge should be correlated. Indeed, the 
monopoles in the Maximal Abelian gauge are by construction Abelian objects, the field of 
which have a natural electromagnetic interpretation. In particular, in order to explain the 
confinement phenomenon as due to these Abelian monopoles, they should be condensed 
in the vacuum of gluodynamics. Contrary to that, the Dirac strings and monopoles, 
Eqs. (^,^), are non-Abelian field configurations. There is no simple electromagnetic- 
like interpretation of them and therefore their relation to the confinement is a separate 
question which should be considered anew. In particular, the contribution of the field 
configurations which correspond the monopoles (|46|) to the expectation value of the Wilson 
loop is unknown. Moreover, one cannot say that if the monopoles ( ^6|) are related to the 
confinement, they should be condensed in the vacuum. Therefore, the correlation of the 
Abelian monopoles in the Maximal Abelian gauge and the monopoles defined by Eq. (^) 
is an interesting question, which we leave for future investigation. 

Due to the non-Abelian nature of the strings (|45|) it is not evident a priori that they are 
unphysical. Indeed, in case of Abelian Dirac strings the U(l) gauge freedom guarantees 
that no physical result depends on the string position. This is true in particular in 
Abelian gauges of gluodynamics. Contrary to that, the strings ( P^D are considered in 
the Lorenz gauge in which there is no remaining U(l) gauge symmetry. Therefore the 



physical irrelevance of the strings (|45|) should be checked separately. 

We have numerically found that, indeed, the strings (^) are unphysical objects. For 
example, their density, defined analogously to the Eq. (|48|) , scales at weak coupling ac- 



cording to the Eq. (PUj), whereas the density of physically relevant string-like excitations 

6n^f} _|_ 102 
11 121 



behaves as Pstring ~ exp{ — + Y^ln/3}. Furthermore, the action density calculated 



13 



around the string away from its boundary is the same as the vacuum action density within 
the numerical errors. Thus the strings ( ^5]) are similar to the Dirac strings in the compact 
electrodynamics. The work is currently in progress to find the rigorous explanation of 
this fact. 



Conclusions. 

We have considered the well known quantum mechanical notion of the geometrical 
Berry phase in the context of non Abelian gauge theories. The bridge between the Berry 
phase and non- Abelian gauge theories is provided by the Wilson loop which can be con- 
sidered as a quantum mechanical evolution operator. The role of the corresponding time- 
dependent Hamiltonian is played then by the gauge potential A{t). This observation 
allows to utilize in gauge theories results obtained earlier for the Berry phase in case of 
Hamiltonians belonging to a Lie algebra of group G. In particular, it is known that the 
space of physically distinct states in this case is the coset space G/H, H being the Cartan 
subgroup of G. Moreover, the evolution of this quantum mechanical system is naturally 
described by the generalized coherent states. 

We have applied these ideas to the Wilson loop operator in the SU(2) gluodynamics 
both in the continuum and lattice formulations. In particular, we have presented an 
explicit construction of the phase ip of the Wilson loop, which is not restricted to the 
interval [— 7r;7r], but still satisfies 1/2TtW = cosy). We have also shown that in the 
coherent states basis the phase ip naturally decomposes into dynamical and Berry phases 
and discussed the topological construction inherent to the Berry phase in gluodynamics. 
Our considerations, when applied to an elementary plaquette, give a lattice analog of the 
Dirac strings and monopoles. The U(l) group which is usually associated with monopole- 
like defects is now defined locally in terms of the infinitesimal Wilson loops, or elementary 
plaquettes. 

The division of the total phase into the dynamical and geometric phases is gauge de- 
pendent. In case of quantum mechanics such a decomposition also depends on a particular 
choice of the basis. We argue, therefore, that for physical applications one has to consider 
a gauge in which the gauge fields are as smooth as possible, thus suppressing all gauge 
dependent singularities. Following the logic similar to that of the paper in Ref. |jll| we 
identify this gauge as the Lorenz gauge. 

We have studied the dynamics of the newly defined monopoles in numerical simulations 
on the lattice and presented a strong evidences that these objects are relevant in the 
continuum limit. Since in the continuum limit the field configurations which produce 
these monopoles become singular, our findings support for the idea that even in the 
Lorenz gauge, in which the gauge potentials are as smooth as possible, the singular fields 
are important. More generally, our results provide with a measure of singular fields in the 
vacuum. Phenomenology related to such fields in the continuum gluodynamics has been 
discussed in Ref. [|l^ 

It is worth emphasizing that the construction we have presented is essentially non- 
Abelian. Only in the particular case when all the gauge fields have the same color ori- 
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entation, it reduces to the known construction of Dirac strings and monopoles in the 
compact U(l) gauge model. Therefore, the relevance of our considerations to the con- 
finement requires further investigation. It is also important to study the relation, if any, 
between monopoles we have constructed and other topological defects, for example, the 
Abelian monopoles in the Maximal Abelian gauge of the gluodynamics. 
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Figure 1: The logarithm of the monopole density, Eq. (^), as a function of /3 = 4/(7^. (a) 
- the symmetric 12^ lattice; (b) - the asymmetric 4 x 12^ lattice. Solid curve represents 
the renormalization group prediction, Eq. (|49D . 
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Figure 2: The asymetry of the 

monopole currents, Eq. versus 
f3 on the 12^ (diamonds) and 4 x 12^ 
(squares) lattices. 



Figure 3: The relative access of 
the magnetic action density, Eq. (0), 
around the monopole currents, Eq. (^), 
for the 12*^ lattice. 



